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The experiments on verification of the Kibble-Zurek mechanism showed that topological defects 
are formed most efficiently in the systems of small size or low (quasi-)dimensionality, whereas in the 
macroscopic two- and three-dimensional samples a concentration of the defects, as a rule, is strongly 
suppressed. A reason for universality of such behavior can be revealed by considering a strongly- 
nonequilibrium symmetry-breaking phase transition in the simplest (p* field model. It is shown that 
the resulting distribution of the defects (domain walls) is formally reduced to the well-known Ising 
model, whose behavior changes dramatically in passing from a finite to infinite size of the system 
and from the low (D=l) to higher (D>2) dimensionality. 
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I. INTRODUCTION 
A. The concept of topological defects 

Formation of topological defects by the strongly- 
nonequilibrium symmetry-breaking phase transitions is 
the subject of interest both in condensed matter and field 
theory. This is because of a close similarity between the 
Lagrangian of Landau-Ginzburg theory, widely used to 
describe phase transitions in the condensed-matter sys- 
tems, and the Lagrangians of the modern elementary- 
particle theories (such as the standard electroweak model 
or various kinds of the Grand Unification Theories), 
which are substantially based on the concept of spon- 
taneous symmetry breaking. 

After the phase transitions in all the above-mentioned 
cases, the stable topological defects of the order param- 
eter can arise, such as the monopoles, strings (vortices), 
and domain walls, depending on the symmetry group in- 
volved. 

The most efficient mechanism of the defect formation 
is the so-called Kibble-Zurek scenario, which is based 
on the simple causality arguments. Namely, if during the 
phase transition the information about the order param- 
eter can spread over the distance ^eff, then the phases of 
the order parameter should be established independently 
in the regions of characteristic size ^cff*^ As a result, af- 
ter some relaxation following the phase transition, the 
stable defects of the order parameter should be formed 
at the typical separation ^cff from each other. So, their 
concentration in a 3-dimensional system can be roughly 
estimated as 



(1) 



where d = 3, 2, and 1 for the monopoles, strings (vor- 
tices), and domain walls, respectively; while the effective 



correlation length foff depends on the particular system 
under consideration. 

For example, in the symmetry breaking of Higgs fields 
by the cosmological phase transitions (i.e. a generation of 
mass of the elementary particles) ^cff is commonly taken 
to be? 



Ccff < c/H, 



pt : 



(2) 



where c is the speed of light, and Hpt is Hubble constant 
at the instant of phase transition. 

In the consideration of vortex generation by a super- 
fiuid phase transition in helium, the corresponding quan- 
tity is 



off 



C2Tq , 



(3) 



where C2 is the speed of the second sound (i.e. a char- 
acteristic speed of propagation of information about the 
phase of the order parameter), and Tq is the so-called 
quench time (i.e. a characteristic time of the phase tran- 
sition), defined as 



1 dT\ 
T~dt) 



(4) 



The similar definitions of ^cff are used also for other 
systems (superconductors, liquid crystals, etc.). 



B. Review of experimental data 

Although Kibble-Zurek scenario was proposed initially 
in the context of cosmological phase transitions in the 
field theories admitting the symmetry breaking (in fact, 
the first idea of such kind was put forward by BogoliuboAJ^ 
almost 40 years ago), much efforts have been undertaken 
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in the last decade to verify this phenomenon by the lab- 
oratory experiments. These works were originated by 
Chuang, et al^ in 1991, and about a dozen of such ex- 
periments were performed by now in various superfluid 
and superconducting systems. They are summarized in 
Table ^ (This table does not show a number of exper- 
iments with liquid crystals, aimed at the detailed study 
of inner structure of the defects without measuring the 
rates of their formation.) 

Analysis of the table reveals a quite interesting ten- 
dency: the topological dejects are formed most efficiently 
in the systems of small size or low ( quasi- ) dimensional- 
ity^ for example, the one-dimensional multi-Josephson- 
junction loops (MJJL)^^ and annular Josephson tunnel 
j unctions as well as in microscopic hot bubbles of '^He 
produced by neutron irradiation.^'^ On the other hand, 
the concentration of the defects in macroscopic systems 
of higher dimensionality, e.g. the two-dimensional su- 
perconductor filmsiSiiiji^ and three-dimensional volume 
samples of ^He (Refs.|13), was found to he considerably 
less than the theoretical predictions. 

In principle, there may be various explanations of the 
above fact (for example, different material parameters of 
the substances listed in Table P). The aim of the present 
work is to propose another explanation, which is based on 
the universal geometric properties of the systems under 
consideration, namely, their size and dimensionality (or 
quasi-dimensionality) . 

II. THE MODEL OF DEFECT FORMATION 
A. Initial assumptions and equations 



the effective correlation length ^eff, whose definition in 
the Kibble-Zurek scenario was already discussed above 
in Sec. II Al The particular value of ^off is not of im- 
portance here, but we shall assume that it is sufhciently 
large in comparison with a characteristic thickness of the 
domain wall ~ 1//^- As a result, the final pattern of vac- 
uum states after the phase transition will look like a 
distribution of spins on the rectangular grid. 

The key assumption of the Kibble-Zurek mechanism 
is that the final symmetry-broken states of the field f in 
two neighboring cells are completely independent of each 
other. Then, the probability of formation of a domain 
wall between them is given evidently by the ratio of the 
number of statistical configurations involving the domain 
wall to the total number of configurations: 

Pkz = 2/4 = 1/2 ; (9) 

and the resulting concentration of the defects (domain 
walls) will be 

riKz = I ^ ' (10) 

where D is the effective dimensionality of the system. 

In fact, such approach is not sufhciently accurate, be- 
cause a nonzero value of the order parameter formed af- 
ter the phase transition represents a coherent state of 
Bose condensate, in which the specific quantum corre- 
lations may occur even at the distances exceeding the 
effective correlation length <^off derived from the "clas- 
sical" arguments. This kind of correlations was clearly 
demonstrated by the MJJL cxperimenlii^. 



To carry out all calculations in analytic form, let us 
consider the simplest iy9^-model of real scalar field (the 
order parameter), whose Lagrangian 

L{v,t) = \ {{d^^f- (V^)^] - \ [^2- (/.VA)]' (5) 

admits the discrete (Z2) symmetry breaking. 

Two stable vacuum states of this field (which, for the 
sake of convenience, will be marked by the oppositely 
directed arrows) are 

= ±V3o = ±Ai/VA ; (6) 

the structure of a domain wall between them, located at 
X — xq, is described as 

<^(a;) = ±(potanh[-^(x-xo)] ; (7) 
and the specific energy concentrated in this wall equals 



Let a domain structure formed after a strongly- 
nonequilibrium phase transition in this model be approx- 
imated by a regular rectangular grid with a cell size about 



B. Overview of MJJL experiment 

For the sake of completeness, let us briefly remind the 
basic idea and result of the multi- Josephson-junction loop 
experiment by Carmi, et al.^^ whose sketch is presented 
in Fig. [2 

A thin quasi-one-dimensional winding strip was en- 
graved at the boundary between two crystalline grains 
of YBa2Cu307-5 high-temperature superconductor film, 
thereby forming a loop of 214 superconductor segments 
separated by the grain-boundary Josephson junctions. 
This system experienced multiple heating-cooling cycles 
in the temperature range 77 K to ~100 K, which covers 
both the critical temperature of superconducting phase 
transition in the segments of the loop (Tc = 90 K) and 
in the junctions between them (Tcj = 83 85 K). 

There is evidently no order parameter in the entire loop 
as long as T > Tc. Next, when the temperature drops be- 
low Tc but remains above Tcj (i.e. T^j < T < Tc), some 
value of the order parameter should be established in 
each segment, as is schematically shown by arrows in the 
right-hand part of Fig. ^ Since these segments are sepa- 
rated by nonconducting Josephson junctions, it seems to 
be reasonable to assume that the phase jumps between 
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TABLE L Summary of experiments on the defect formation by strongly-nonequilibrium symmetry-breaking phase transitions. 
A positive result (+) implies an agreement with the Kibble-Zurek estimate within a factor about unity; and the negative 
result ( — ), a disagreement by a few times to a few orders of magnitude. 
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FIG. 1: Sketch of the experimental setup (left) and basic results (right) of MJJL experiment. 



them are random (i.e. uncorrelated to each other). Fi- 
nally, when the temperature drops below Tcj, the entire 
loop becomes superconducting and, due to the above- 
mentioned jumps, a phase integral along the loop, in gen- 
eral, should be nonzero. As a result, the electric current 
/, circulating along the loop, and the corresponding mag- 
netic flux penetrating the loop, will be spontaneously 
generated. 

So, if the phase jumps in the intermediate state Tcj < 
T < Tc were absolutely uncorrelated, then distribution of 
the spontaneously-trapped magnetic flux in the particu- 
lar experimental setupM. would be given by the normal 
(Gaussian) law with a characteristic width 3.6 (j)o (where 
(j)o is the magnetic-flux quantum). On the other hand, 
the actual experimental distribution was found to be over 
two times wider; and this anomaly was satisfactorily in- 
terpreted by the authors of the experiment assuming that 



the phase jumps in the intermediate state were not ran- 
dom but correlated to each other so that the probability 
P{Si) of the phase difference 5i in the i'th junction was 



Pi6i)^exp[~Ej{d,)/T,] 



(11) 



where Ej is the energy concentrated in the Josephson 
junction, and is the phase transition temperature, 
measured in energy units. 

The nature of the above-mentioned correlations be- 
tween the phases of order parameter in the spatially- 
separated regions is not sufficiently understood yet. They 
may be one of manifestations of the specific quantum "en- 
tanglement" , revealed by now in various quantum sys- 
tems. We are not going to discuss here this problem in 
more detail but assume that the correlations like Hll|) are 
typical for all Bose condensates formed by the symmetry- 
breaking phase transitions. As will be seen from the sub- 
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TABLE II: Formal correspondence between the phase tran- 
sitions in ifi* and Ising models. 



Ising model 



1. State formed after a 
strongly-nonequilibrium 
phase transition. 

2. Elementary domains of 
the symmetry-broken state. 

3. Statistical sum Z for 
distribution of the domains 
with various vacuum states. 

4. Critical temperatures Tc 
for various systems. 

5. Suppression of the domain 
wall formation. 



Thermodynamically- 
equilibrium state. 

Spins. 

Thermodynamical statistical 
sum Z for spin distribution. 

Temperature T varying in 
the course of the phase 
transition. 

Phase transition in the spin 
system. 



sequent analysis, such an assumption will enable us to 
interpret satisfactorily a wide range of experimental data 
on the efficiency of formation of topological defects. Just 
this fact can be considered as indirect confirmation of 
universality of Hll|l . 



C. Improvement of the standard estimates 

As follows from the above discussion, the probabilities 
of various field configurations after a symmetry-breaking 
phase transition should be calculated taking into ac- 
count the energy concentrated in the defects. As a re- 
sult, a pattern of the domains formed by the strongly- 
nonequilibrium phase transition in the tp^ lattice model 
will look like a distribution of spins in Ising model at the 
temperature T which is formally equivalent to the criti- 
cal temperature Tc of the initial (^^-model. All aspects of 
this formal correspondence are summarized in Table UTI 

Then, the probability of formation of a domain wall 
should be calculated as 



P = 



d 



InZ 



(D) 



(12) 



where E is the energy of the elementary domain wall (i.e. 
at one boundary between two neighboring cells), N is the 
number of cells along each side of the lattice, and 



=^exp(-£,/r,) 



(13) 



is the usual statistical sum over all possible spin config- 
urations of the Ising model, where £i is the total energy 
of j'th configuration. 

Particularly, for 1-dimcnsional system with periodic 
boundary conditions, 

z'^'-Y: E -p|-;^e|(i-^^^^+i)}; (14) 

t=ls,=±l ^ k=l ^ ^ 



for 2-dimensional system, 

N N 

^'^^-EE E -p 



E 
T'c 



N N 

T.T.I 



i = lj = lSij=±l ^ ^ k=l 1 = 1 

SkiSk+1,1 - SkiSk,i+i) I ; (15) 

and so on. Here, Sk and Ski are the spin-like variables 
describing the broken-symmetry states in the fc'th and 
(fcZ)'th cell, respectively. 

As is known (e.g. Ref. Ush . the Ising model for one- 
dimensional as well as the finite-size higher- dimensional 
systems does not experience a phase transition to the or- 
dered state at any value of the ratio E/T. From the 
viewpoint of the domain wall formation by the strongly- 
nonequilibrium phase transition in ip - model (Table mil, 
this means that concentration of the defects will not dif- 
fer considerably from the standard Kibble-Zurek estimate, 
because the probability of defect formation P at the scale 
of the effective correlation length ^cff will not deviate sub- 
stantially from Pkz — 1/2 ~ 1 . 

On the other hand, the Ising model for the sufficiently 
large (infinite-size) two- and three-dimensional systems 
does experience a phase transition to the ordered state 
at some value of E/T ^ 1. As a result, the concentra- 
tion of domain walls in the corresponding ip'^ -model at 
large ratios E/Tc should be suppressed dramatically, due 
to formation of macroscopic regions with the same value 
of the order parameter, covering a great number of cells of 
the effective correlation length ^cff- (To avoid misunder- 
standing, let us emphasize once again that the different 
values of Tc should be understood here as the critical tem- 
peratures for various physical systems described by the 
(yj^-model, and they formally correspond to the variable 
temperature T of a fixed Ising system.) 



D. Particular example 

The general conclusions formulated above can be illus- 
trated by the particular example in Fig. |21 which repre- 
sents the concentration of domain walls n normalized to 
the standard Kibble-Zurek value as function of E /Tc 
for the following three systems: 

(1) one-dimensional Ising model, which admits the exact 
solution (see, for example, Ref. (16 for more details); 

(2) two-dimensional 6x6-cell Ising model with periodic 
boundary conditions, which simulates quite well the sys- 
tem of infinite size; and 

(3) two-dimensional 6x6-cell Ising model with free 
boundaries (where no energy is concentrated), which is 
an example of a microscopic system. 

(The particular size of 6 cells along each side of the lattice 
was taken quite arbitrarily, just as the value at which a 
numerical computation of the statistical sum l|15|l is not 
too cumbersome; details of the calculations can be found 
in Appendix.) 
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TABLE III: Nonzero statistical weight factors for 6x6-cell 
Ising model with periodic boundary conditions. 
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FIG. 2: Concentration of the defects n normalized to the 
standard Kibble-Zurek value Wkz as function of the ratio of 
the domain wall energy E to the phase transition tempera- 
ture Tc for the infinite-size ID Ising model (dashed curve), 
6x6-cell 2D Ising model with periodic boundary conditions 
(solid curve) , and the same model with free boundaries (dot- 
ted curve). 



As is seen in Fig. |21 at E/Tc ~ 1 the concentration 
of defects in one-dimensional system (dashed curve) dif- 
fers from the standard Kibble-Zurek value by less than 
two times; in microscopic two-dimensional system (dot- 
ted curve), by three times; while in the macroscopic two- 
dimensional system (solid curve) it is suppressed by an 
order of magnitude. Such suppression becomes much 
stronger when the ratio E/Tc increases: for example, at 
E/Tc ~ 2 the difference between each of the curves is 
over an order of magnitude. 



III. DISCUSSION 

As follows from the above consideration, the specific 
thermal correlations between the phases of Bose conden- 
sates in separated spatial regions, revealed for the first 
time in MJJL experiment, can show us a general way 
for explanation of the entire diversity of experimental 
data on the topological defect formation. Unfortunately, 
the results of our simplified model cannot be confronted 
quantitatively with the ones cited in Tabled because the 
superfluids and superconductors possess more complex 
order parameters and Lagrangians than So, more 
elaborate calculations are required. 

Next, it is important to emphasize that the refined 
concentration of defects 



f[E/Tc) nKz(Coff<TQ)) 



(16) 



possesses exactly the same dependence on the quench 
rate Tq as in the classical Kibble-Zurek scenario. So, this 
dependence, often measured in the experiments, in gen- 
eral, cannot serve for discrimination between the models; 
the absolute values of the defect concentration are always 
necessary. 



Finally, let us mention that the ideas described in the 
present work can be applied also to solving the problem 
of excessive concentration of topological defects predicted 
after the cosmological phase transitions of Higgs fieldsiii 
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APPENDIX: COMPUTATION OF THE 
STATISTICAL SUMS 

For the sake of completeness, let us present the formu- 
las used for drawing the curves in Fig. |21 

As is known, all characteristics of the one- dimensional 
Ising model can be calculated analytically (for a detailed 
description of the method see, for example, R ef . IT^ . The 
resulting formula for a defect concentration in a chain of 
length N with periodic boundary conditions is 

n/TiKz = 2 exp(-£;/Tc) x (A.l) 

(1 + exp(-i^/r,))^-i - (1 - exp(-g/r,))^-i 
(1 + eM-E/Tc)r + (1 - eM-E/Tc)r ' 

which in the case of infinite chain is reduced to 

2 exp(-£;/T,) 



1 + eM~E/Tc) 



(A.2) 
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TABLE IV: Nonzero statistical weight factors for 6x6-cell 
Ising model with free boundaries. 
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For the two-dimensional infinite lattice, statistical 
sum l)13|l can be reduced by the well-known technique 
to some integral over two variables, and the asymptotics 
of this integral can be obtained analytically just in the 
point of the phase transition of Ising model (for more 
details see, for example, Ref . fisj) . Unfortunately, these 
results are of little value for our work, since we need to 
know the statistical sum in a wide range of temperature. 



Besides, this method cannot be applied to the finite-size 
lattices, which are also of considerable interest for us. 

So, we used a more straightforward approach. In gen- 
eral, the statistical sum (I13|) can be rewritten as 

=Y,Ck eM-kE/ T,) , (A.3) 

k 

where Ck is the statistical weight of configurations in- 
volving k elementary domain walls with energy E each. 
In the particular case of two-dimensional 6x6-cell Ising 
model, the coefficients Ck were calculated by the meth- 
ods of computer algebra both for the lattice with periodic 
boundary conditions (in which every domain wall carries 
the energy E) and with the free boundaries (where only 
the inner domain walls contribute to the total energy of 
the configuration). The corresponding coefficients C^°' 
and C^''°<= are listed in Tables UHl and Hvl 

Next, the concentration of domain walls is obtained by 
substituting (|A.3|I into H12|) ; and the final result takes the 
form 

2 V kCkexpi-kE/Tc) 

n/"Kz = , (A.4) 

D7V°^^Cfcexp(-fc£;/r,) 

where in the case under consideration D = 2 and 6. 
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